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Abstract. I review some key aspects of the ongoing progress in our understanding of the 
infrared dynamics of the QCD Green’s functions, derived from the close synergy between 
Schwinger-Dyson equations and lattice simulations. Particular attention is dedicated to the 
elaborate nonperturbative mechanisms that endow the fundamental degrees of freedom (quarks 
and gluons) with dynamical masses. In addition, the recently established connection between 
the effective interaction obtained from the gauge sector of the theory and that needed for the 
veracious description of the ground-state properties of hadrons is briefly presented. 


1. Introduction 

Sixty years after their invention [T] , Yang-Mills theories occupy the center stage of elementary 
particle physics, providing a fundamental description for an impressive array of physical 
phenomena. The remarkable property of asymptotic freedom makes these theories particularly 
attractive, and, in a sense, self-contained. However, while perturbation theory works well in 
their ultraviolet regime, the infrared dynamics represent a challenge of notorious complexity. 
In fact, in the case of QCD [2], the most characteristic phenomena, such as confinement, mass 
generation, and bound-state formation, are purely nonperturbative. 

In this presentation we will focus on the paradigm-shifting picture that emerges for some of 
the aforementioned phenomena from the systematic study of Green’s functions within a global 
framework, where the analysis of the Schwinger-Dyson equations (SDEs) is complemented and 
refined by inputs obtained from large-volume lattice simulations. 

The present study is restricted to the special case of the Landau gauge, where the vast 
majority of recent lattice simulations have been performed. The SDE part of the problem 
is addressed within the formalism obtained from the fusion between the pinch technique 
(PT) [SHUEIEIIT] and the background field method (BFM) [8], denoted simply as “PT-BFM” (for 
a representative sample of different approaches, see p iini ini iia [la [la [la [ig iia iig iig [2^ ). 

This presentation is organized as follows. In section[2]we review the PT-BFM formalism, and 
highlight its most salient features. Then, in section [3] we summarize the subtle filed-theoretic 
mechanism that leads to the gauge-invariant generation of a dynamical gluon mass. In section [3] 
we present a brief overview of the impressive coincidence achieved between SDE calculations and 
lattice simulations of the ghost propagator, while in section [5] we address some of the intricate 
issues surrounding the dynamical generation of the constituent quark masses. Then, in section [6] 
we review a recent ab-initio derivation of the interaction kernel that has been extensively used 



Figure 1. The conventional SDE of the standard gluon propagator (QQ). Black blobs 
represent fully dressed one-particle irreducible vertices, whereas the white ones denote fully 
dressed propagators. 


in the contemporary description of hadron phenomenology. Finally, in section [7] we summarize 
our conclusions. 

2. The PT-BFM formalism: a powerful framework for studying SDEs 

In what follows we will work exclusively in the Landau gauge, where the gluon propagator 
assumes the totally transverse form 

)) ~ 9^iv ~Qij.Qu/q ■ ( 1 ) 

The scalar form factor A(g^) is related to the all-order gluon self-energy n^i,(g). Specifically, 
as a consequence of the BRST symmetry, is transverse, both perturbatively and 

nonperturbatively; one has then 

q^U^uiq) = 0 ; ( 2 ) 

and = q‘^ + 17(5^). In addition, for later convenience, we introduce the inverse of the 

gluon dressing function, to be denoted by J{q'^), namely 

A-\q^) = q^J{q‘^). (3) 

The nonperturbative dynamics of the gluon propagator are governed by the corresponding 
SDE. In particular, within the conventional formulation [ZI1I22], the gluon self-energy Il^,^{q) 
is given by the fully dressed diagrams shown in Fig. [TJ One longstanding difficulty with this 
equation is that it cannot be truncated in any obvious way without compromising the validity of 
Eq. ([2|). This happens because the conventional fully dressed vertices appearing in the diagrams 
of Fig. [1] satisfy complicated Slavnov-Taylor identities (STIs), and it is only after the inclusion 
of all diagrams that Eq. ([2]) may be enforced. 

Instead, the formulation of this SDE in the context of the PT-BFM formalism furnishes 
considerable advantages, because it allows for a systematic truncation that respects manifestly, 
and at every step, the crucial identity of Eq. ([2]). The deeper reason why such a truncation 
becomes possible may be traced back to the fact that the PT-BFM Green’s functions satisfy (by 
construction) Abelian-like Ward identities (WIs) instead of STIs. In particular, if we employ 
the BFM language and distinguish the gluon helds into background {B) and quantum (Q) 



Figure 2. The SDE obeyed by the QB gluon propagator. Black blobs represents fully dressed 
1-PI vertices; the small gray circles appearing on the external legs (entering from the right, only!) 
are used to indicate background gluons. The diagrams contained in each box form individually 
transverse subsets. 


ones, vertices such as BQQ, Bcc, and BQQQ, to be denoted by Ta^u{Q,'f',p) , Ta{q,r, —p), and 
T'^i^p^{q,r,p,t), respectively, satisfy the simple WIs [23] 

{q, r, p) = il^~l (r) - iA'J (p), (4) 

(?“f r, -p) = D-\q + r) - D-\r), (5) 

and 

= n^r'^T,p^{r,p,q + t) + r^^r^T,,,{p,t,q + r) 

+ n^nra.p{t,r,q+p). ( 6 ) 

These particular vertices (instead of the conventional ones) appear in the SDE of Eig. [2l 
which controls the self-energy Il^i,(q) of the the mixed background-quantum gluon propagator 
(BQ), denoted by A(g^). Consequently, it is relatively straightforward to prove the block-wise 
transversality of n^iy(q), namely [53] 

9^[(«i) + (« 2 )]mi/ = 0; q'^iias) + ( 04 )]^!^ = 0; ^^[( 05 ) -h (og)];,!. = 0. (7) 

Erom Eq. ([2|) it is clear that the SDE of A(g^) contains inside its defining diagrams the full 
propagator A(g^); therefore, in that sense, it cannot be considered as a dynamical equation for 
neither A(g^) nor A(g^). At this point a crucial identity relating A(g^) and A(g^) |241I25| enters 
into the game. Specifically, one has that 

A(g2) = [l + G(g2)]A(g2), (8) 

where G{q^) is the co-factor of a certain two-point function that originates from the ghost 
sector of the theory. 

The fundamental observation put forth in a series of works [231123 Cn is that one may use 
the SDE for A(g^), take advantage of its improved truncation properties, and then convert it 
to an equivalent equation for A(g^) (the propagator simulated on the lattice) by means of the 









special relation given in Eq. ([8]). Thus, the PT-BFM version of the SDE for the conventional 
gluon propagator reads 




1 + G(g2) 


(9) 


3. Dynamical gluon mass generation 

The possibility that Yang-Mills theories generate dynamically a gluon mass was first proposed 
and explored in the early eighties O EHl EH], and has attracted particular attention in 
recent years, mainly due an accumulation of hard evidence obtained from large-volume lattice 
simulations, both in SU(3) (SO] [3ll [32l [33] and in SU(2) [3H [35l (36] [37] (for a different lattice 
approach, see [55]h 

What these simulations reveal is that the Landau gauge gluon propagator saturates in the 
deep infrared at a fixed nonvanishing value, as shown on the left panel of Fig. [S] which is the 
smoking gun signal for gluon mass generation. A conclusive demonstration of how this may 
happen at the level of the gluon SDE within the PT-BFM framework has been given in |39| . 



Figure 3. The quenched lattice gluon propagator A{q‘^) (left panel), and the corresponding 
gluon mass (right panel). 

From the kinematic point of view we will describe the transition from a massless to a massive 
gluon propagator by implementing the replacement (in Euclidean space) 

A-\q‘^)=q^J{q‘^) A-\q^) = q^J{q'^) + m\q^), ( 10 ) 

where m?{q'^) is the (momentum-dependent) dynamically generated mass, with the property 
that m^(0) > 0; evidently, A”^(0) = m^(0). 

The actual held-theoretic mechanism that permits the generation of the term m^(g^) can 
be traced back to the seminal work of Schwinger [Ml SI], which may be summarized in the 
statement that the vacuum polarization of a gauge boson that is massless at the level of the 
original Lagrangian may develop a massless pole, whose residue will be eventually identified 
with m^(0). In the case of Yang-Mills theories, the origin of the aforementioned poles is due to 
purely nonperturbative dynamics: for sufficiently strong binding, the mass of certain (colored) 
bound states may be reduced to zero [izisaiiiisiiB]. As has been shown in m through a 
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Figure 4. The SDEs for the full three and four gluon vertices in the presence of their 
corresponding pole parts. 


detailed study of the Bethe-Salpeter equation that controls the formation of such a composite 
massless pole, this dynamical scenario is indeed realized. 

In addition to triggering the Schwinger mechanism, these bound-state poles act as composite, 
longitudinally coupled Nambu-Goldstone bosons, maintaining gauge invariance; notice, however, 
that they differ from ordinary Nambu-Goldstone bosons as far as their origin is concerned, 
since they are not associated with the spontaneous breaking of any continuous symmetry. In 
particular, gauge invariance requires that the replacement described schematically in Eq. m 
be accompanied by a simultaneous replacement of all relevant vertices by 

r —^ f' = f + E, (11) 

where V contains the aforementioned massless poles (see Fig. 11, and guarantees that the new 
vertex E' satisfies the same WIs as E, but now replacing the gluon propagators appearing on 
their rhs by massive ones. For example, in order to maintain Eq. (H intact, the corresponding 
V must satisfy the WI 


q'^Vaf,uiq,r,p) = - rr?{p^)P^,u{p)] ( 12 ) 

similarly, when contracted with respect to the other two momenta, Vafiu{q, r,p) satisfies the STIs 
necessary to enforce the validity of the corresponding STIs satisfied by Ta^^{q,r,p), but now 
with massive instead of massless gluon propagators on the rhs. 

The fact that the massless poles must be longitudinally coupled is reflected at the level of T' 
through the condition 

P^'^iq)P^^'^{r)P’''’^{p)V^,^,,,{q,r,p) = 0, (13) 

and a exactly analogous relation for the four-gluon Va^upiq,r,p,t). These latter conditions 
guarantee that the massless poles deconple completely from the on-shell 5-matrix or other 
physical observables. Moreover, when Eq. (jl3p is combined with Eq. m and the additional 
STIs not reported here, one obtains the exact closed expression for Va'p'u'{q,'f',p) [i8] . 


The inclusion of the aforementioned vertices into the gluon SDE leads to its separation into 
two integral equations of the generic form 

Jm{q^) = 1 + / 

Jk 

= [ IC2{q‘^,m‘^,Am), (14) 

Jk 

such that q'^}Ci{q‘^, m?, Am) —>■ 0, as 0, whereas }C 2 {q‘^, m?, Am) / 0 in the same limit, 

precisely because it includes the 1/g^ terms contained inside the V terms. This set of equations 
bares a close analogy to those studied in the case of the more familiar phenomenon of quark 
mass generation; in particular, m^{(f‘) plays the role of the quark mass function (denoted by 
B{p^) in Eq. (f^ b 

The explicit closed form of the integral equation that determines m?‘{(f‘), together with a 
variety of related theoretical considerations, have been presented in a series of works gaEolEI]. 
The upshot is that one obtains a dynamical gluon mass of the form shown on the right panel 
of Eig. [21 whose functional dependence may be fitted by 


^ 1 + ((?2/A42)i+P 


(15) 


with mo = 375 MeV, p = 0.11 and M. = 431 MeV. 

The nontrivial momentum dependence of the gluon mass is mainly responsible for the fact 
that, contrary to a propagator with a constant mass, the gluon propagator of Eig. [3] displays an 
inflection point. It turns out that the presence of such a feature is a sufficient condition for the 
spectral density of the gluon propagator, p, to be non-positive definite. 

Specifically, the Kallen-Lehman representation of the gluon propagator reads 

and if A(g2) has an inflection point at g*, then its second derivative vanishes at that point 
(see Fig. B, namely [52] 

A"(g^) = 2^ =0- (17) 



Figure 5. The first and second derivatives of the gluon propagator. 










Given that > 0, then p{a) is forced to reverse sign at least once. The non-positivity of p{cr), 
in turn, may be interpreted as a signal of confinement (see [53], and references therein), since 
the breeching of the axiom of reflection positivity excludes the gluon from the Hilbert space of 
observables states (for related works, see [Ml ISSl [56l [571 ISEl ISHl [60] )• 

As can be seen in Fig. [5] the first derivative of A(g^) displays a minimum at q‘^ = 0.238 GeV^, 
and, consequently, the second derivative vanishes at that same point. 

4. Ghost sector 

The nonperturbative behavior of the ghost propagator, D{p^), has been the focal point of 
intensive studies, because in the past it had been linked to a certain confinement scenario [9], 
whose realization required that the corresponding ghost dressing function, F{p^) = p^D{p^), 
should diverge at = 0. However, a plethora of high quality lattice simulations (for 
SU(3) [211 EH EQI [33] and for SU(2) [HU [35l [36l EZ] ) have conclusively established that F{p‘^) is 
finite at the origin. In fact, the finiteness of F{p^) is tightly interwoven with the massiveness 
of the gluon, because it is precisely the presence of the gluon mass that tames the logarithm 
associated with T(p^), and prevents it from diverging in the infrared [39]. 



Figure 6. The SDE for the ghost propagator given by Eq. (|18p . The white blobs represent the 
fully dressed gluon and ghost propagators, while the black blob denotes the dressed ghost-gluon 
vertex. 

Even though the finiteness of Fijp') may be explained qualitatively in the above relatively 
simple terms, a more strenuous effort is required in order to obtain from the corresponding SDE 
the full curve of F{p^) observed on the lattice. 

The SDE for the ghost propagator is diagrammatically represented in the Eig. [B] and assumes 
the form 


W \p^) =ip^ - q^Ca 



,-k - p,p)l^^''{k)Ty{-k, -p, k+p)D{k+p). (18) 


In the above equation Ca represents the Casimir eigenvalue of the adjoint representation {N for 
SU{N)), d = 4 — e is the space-time dimension, and we have introduced the integral measure 




(27r)« 


dU, 


(19) 


with p the’t Hooft mass. The vertex Ti^{—k, —p, k + p) is the fully dressed ghost-gluon vertex, 
whose tensorial decomposition is given by {r = k + p) 


T^{-k, -p, r) = A{-k, -p, r) py + B{-k, -p, r) ky ; 


( 20 ) 











Figure 7. Left panel: The numerical result for the form factor A{—k, 0, k) in the soft ghost 
configuration (red continuous line) [6l!- Right panel: The F{p^) obtained from the ghost SDE 
with a bare ghost-gluon vertex (blue dashed curve) and the corresponding solution when the 
A{—k,0,k) is used (red continuous line) , using in both cases jjL = 4.3 GeV and as{p) = 0.22. 
The lattice data are from Ref. m- 


at tree-level, the two form factors assume the values A^^\—k, —p, r) = 1 and B^^\—k, —p, r) = 0, 
reproducing the standard bare vertex = p^. Clearly, due to the full transversality of A^j^(fc), 
any reference to the form factor B disappears from the ghost SDE of Eq. (1181) . and one obtains. 


F-\p^) 


1 -|- ig'^CA [ 
Jk 



jk-p) 

J^2p2 


21 


A{—k, —p, k + p)A{k)D{k + p). 


( 21 ) 


The dependence of this equation on both A(/c) and A{—k, —p,k + p) is just another example 
of the perennial difficulty inherent to the SDE framework, namely the fact that, in principle, 
all Green’s functions of the theory are coupled to each other. In the case at hand, A{k) is 
known from the lattice (Fig. [3]), and there is no need to resort to its own SDE. However, A is 
more problematic, since there are no lattice data for the kinematic configuration relevant to this 
problem (for SU(2), see [UTl [S2l E3])- The obvious easy way out, namely the use of the tree- 
level value for A is quantitatively insufficient: one obtains the result shown on the right panel 
of Fig. [7] (blue dashed curve), which, even though qualitatively correct, it is clearly considerably 
suppressed compared to the lattice data. 

To ameliorate this situation, an approximate version of the SDE equation that controls the 
evolution of A has been considered in [M], in the special kinematic limit p —)• 0 (for a related 
study, see my, the result for A{—k, 0, k) is shown in the left panel of Eig. [71 

Even though the increase with respect to the tree level value is not too dramatic, the peak 
of A{—k,0,k) at around 1 GeV turns out to be just right: when coupled with Eq. (|2T]) . a 
considerable enhancement of the resulting Fljp') is produced, due to the nonlinear nature of this 
system of integral equations. As shown in the right panel of Eig. [71 one may actually reproduce 
the lattice data very accurately, using the correct value for the strong coupling Os [Ml EZ] • 


5. Quark sector 

It is well known that the nonperturbative mechanism that endows quarks with their 
constituent masses accounts for about 98% of the visible matter in the Universe. In fact, it 










has often been denominated as the most efficient mass generating mechanism known, practically 
generating masses from nothing (see [53] and references therein). This fundamental phenomenon 
(and the related dynamical breaking of the chiral symmetry) is encoded in the SDE of the quark 
propagator, know also as quark “gap equation”. 

The standard decomposition of the quark propagator, S{p), is given by [2T] 

S~^{p) = A{p‘^)^ - B{p'^)I = A{p‘^)\^ - M{p‘^)l], (22) 

where I is the identity matrix, and }A{p^) is the dynamical quark mass; its generation is 
inextricably connected with the dynamical breaking of the chiral symmetry. 


^ — TTlo + (->-) ^ + 

Figure 8. The SDE (gap equation) for the quark propagator, given by Eq. (|23l) . The white 
blobs represent the fully dressed gluon and quark propagators, while the black blob denotes the 
dressed quark-gluon vertex. 

The SDE for the quark propagator is shown diagrammatically in Fig. (8] and is given by 

S-\p) ^ rp5(A:)r,(-p, k, q)A^^{q) , (23) 

where q = p — k, r,y{—p, k, q) is the fully dressed quark-gluon vertex, is the tree-level value, 
and Cp is the Casimir eigenvalue in the fundamental representation. In this talk we will only 
consider the case of vanishing current quark mass, niQ = 0, i.e., the chiral symmetry is kept 
intact at the Lagrangian level. One of the characteristic features of this integral equation is that, 
in order to give rise to nontrivial solutions, the support of its kernel in a special momentum 
regime must overcome a critical amount. 

As happens in the case of the ghost equation, the nonperturbative vertex constitutes a 
crucial ingredient for overcoming the aforementioned critical amount and arriving to a physically 
acceptable solution. The situation, however, is much more complicated now, given that the 
tensorial decomposition oiT^{—p,k,q) consists of twelve independent structures [72] instead of 
only two. Projecting them out of the SDE that governs Ti,{—p, k, q) represents a major technical 
challenge, let alone solving the resulting intricate system of coupled integral equations. Lattice 
simulations, on the other had, have obtained some of the corresponding form factors for special 
kinematic configurations but their relevance to the problem at hand is relatively 

limited. 

The alternative that has been extensively adopted is to employ gauge-technique inspired 
Ansatze (see, e.g., m and references therein) for Ti,(—p,k,q), namely construct it in such a 
way that it satisfies the STI imposed by the BRST symmetry. In particular, T^{pi,p 2 ,P 3 ) obeys 
the fundamental STI |2| 

P3'^ti{Pl,P2,P3) = F{P3)[S~\-PI)H{PI,P2,P3) - H{P2,P1,P3)S~^{P2)] , (24) 

where the fermion-ghost scattering kernel H(pi,p 2 ,P 3 ) is defined diagrammatically in Fig. (9] 
and has the tensorial decomposition 



H{pi,P2,P3) = +X2:^2 + X3a^yp'^pl, 


(25) 







Figure 9. Diagrammatic representation of the ghost-gluon kernel H. 


with H{p 2 ,Pi,P 3 ) is the “conjugated” version of H{pi,p 2 ,P 3 ) UM- 

Note that when both F and H are set to unity (i.e., the ghost sector is switched off), Eq. (I2^ 
reduces to the text-book WI relating the photon-electron vertex with the electron propagator 
in QED. The the so-called Ball-Chiu (BC) vertex [73], 




A{pi) + A{p2) 


-T 


+ 


{pi -P 2 Y 

p\-pI 


[A{pi) - A{p2)\ 


+ [B{pi) - B{p2)] 


(26) 


satisfies precisely this particular WI, instead of the full STI of Eq. and has served as the 
starting point of numerous studies. 

Specifically, use of Eq. ([261) into Eq. (1231) leads to a coupled system for A{p‘^) and B{p^) 


A{pY 1 + Cpg'^ ^ A^k^k^ + ’ 

B{pY = A2(fc2)l^ + B2(A:2)^^^^^’^^’ 

where the kernel }C^^{k,p) originates from Eq. (|26p (for its explicit form, see, e.g., |74]i. while 
the kernel JC{q) encompasses all remaining contributions. 

If the ghost contributions to E^ are switched off, JC{q) consists simply of A(g); in this case, 
however, the resulting combined kernel (right panel of Fig. fTOll fails to furnish a nontrivial 
solution (left panel of Fig. fTOll . 

Given this unsatisfactory situation, a next possible improvement consists in restoring the 
omitted ghost-related contributions to E^, thus achieving the gradual “non-Abelianization” of 
the Ball-Chiu vertex. In particular, one may include the contributions of F(q) and F[, as dictated 
by the STI of Eq. (|24ll . Evidently, since F(q) is very well known (see previous section), the main 
challenge is to obtain an approximation for F[, or at least, for some of its form factors. In 
addition, renormalization introduces (effectively) into JC{q) an additional multiplicative factor 
of F{q), an approximation which guarantees that the anomalous dimension of the resulting 
dynamical quark mass is consistent with the known one-loop expression (see also m)- On the 
right panel of Fig. [TOlone may follow the gradual built-up of IC{q) and on the left the form of the 
quark mass M.{p^) obtained at each step. Evidently, when K,{q) = F‘^{q)A{q), a considerable 
amount of M.{p^) is obtained from Eq. (l271l : however, one deviates by a factor of about 30% 
from the phenomenologically accepted value. 

It turns out that the missing amount of dynamical quark mass may be obtained by supplying 
to /C(g) some of the strength coming from H. In particular, the form factor Xq was computed 
within the one-loop dressed approximation, for the symmetric kinematic configuration 

{—q/2, —ql2, q) [71|. The result is shown in Fig. [TTl note that the deviation from unity (peaked 
around 600 MeV) is modest, but, due to the nonlinearity of the gap equation, it gives rise to a 
considerable amount of additional quark mass, as can be appreciated in Fig. [TOl 








Figure 10. Left panel: The individual contribution of the ingredients composing K,{q). 
The green dotted-dashed line represents the case where lC{q) = A(g), the blue dotted line 
the case where K,{q) = A{q)F{q), in the red dashed curve lC{q) = A{q)F‘^{q) and, finally the 
black continuous line represents the case where }C{q) assumes the full form used, namely 
}C{q) = A{q)F‘^{q)X^\q). Right panel: The corresponding dynamical quark mass generated 
when we use in Eq. (1271) the different forms of JC{q) shown in the left panel. 



Figure 11. Numerical result for the form factor Xg^^(( 7 ), obtained within the one-loop dressed 
approximation, using fully dressed propagators and bare vertices m- 


6. Hadron observables from continuum QCD 

One of the longstanding challenges of QCD is to furnish quantitatively accurate ab initio 
predictions for the observable properties of hadrons [761 EZl ESI EH]. In a recent work |8n| . 
considerable progress has been made in this direction; in what follows we will highlight the main 
aspects of these new developments. 

Let us consider the kernel appearing in a typical Bethe-Salpeter equation [8T], shown 
in Fig. fT2l which is contained in the blue box. This kernel receives a “universal” (process- 
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Figure 12. A typical Bethe-Salpeter equation. The blue box contains the interaction kernel, 
whose form is crucial for the type of solutions (spectrum) obtained. 


independent) contribution, whose origin is the pure gauge sector of the theory. In that sense, 
this contribution constitutes the common ingredient of any such kernel, regardless of the nature 
of the particles between which it is embedded; most notably, it does not depend on the valence- 
quark content of the Bethe-Salpeter equation. 

The systematic diagrammatic identification of the precise pieces that constitute this particular 
quantity may be carried out following the PT. As has been explained in detail in the literature, 
the upshot of this construction is the rearrangement of a physical amplitude into sub-amplitudes 
with very special properties [?]• In addition, it is well-known that the result of this particular 
rearrangement coincides with the BB gluon propagator dehned in the BFM (see [6] and 
references therein). 

Specifically, the standard A(q^) introduced in Eq. ([1]) and the corresponding quantity of the 
BB gluon propagator, denoted by A(q^), are related by the exact relation 

A(g2) = A(g2)[i + G'(g2)]2, (28) 


which is completely analogous to Eq. (jSj). 

At the one-loop level, and keeping only ultraviolet logarithms, one has |39| 


1 + G(g2) 


A-\q^) 



1 + 


oisCa 

Utt 


In 




oisCa 

Utt 
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(29) 


and thus 


A-\q^) = q^ 


1 -I- bag 



(30) 


where b = IICa/IStt is the first coefficient of the Yang-Mills /3 function, as it should [8]. 

Due to the special WIs satisfied by the PT-BFM Green’s functions, the (dimensionful) 
universal combination 


d{q^) = asA{q‘^) 


asAjq^) 

[1 + G(rf 


(31) 


is renormalization-group invariant. In particular, the renormalization constants of the gauge¬ 
coupling and of A~^(p^), dehned as 


9{P^) = Zg 

Aip^g^) = Zj\p^)Ao{p^), 


( 32 ) 





















Figure 13. Left panel:: The gluon propagator for different values of the renormalization point 
pi. Right panel: The auxiliary function 1 + G{q^) for the same set of values of p. 


where the “0” subscript indicates bare quantities, satisfy the QED-like relation 

(33) 

Therefore 

do{q^) = asoAo(g^) = as{fi‘^)A{q‘^, = d{q^) (34) 

maintains the same form before and after renormalization, i.e., it forms a /r-independent 
quantity. 

This last property has been explicitly verified to an impressive accuracy, following the 
procedure summarized in Fig. [T3] and Fig. fT41 The gluon propagator A(g^), obtained from 
the lattice data of [31], is multiplicatively renormalized (in the momentum-subtraction scheme) 
at hve different values of p. The dynamical equation furnishing 1 -|- G{q^) (not reported here, 
see [82l[83|) is subsequently solved, for the same values of p. When the resulting curves (given 
in the two panels of Fig. [T3ll are combined according to Fq. m, it turns out that a unique 
curve emerges, provided that the values of each as{p^) are those shown in the insert of Fig. [TTl 
From the dimensionful quantity d{q^) one may extract a dimensionless effective charge, to be 
denoted by T(g^), by simply pulling out a factor of q^, namely 

Z{q^) = q^d{q^). (35) 


This SDF-derived quantity allows for a direct comparison with an analogous quantity 
commonly employed by practitioners for describing the (momentum-dependent) interaction 
between quarks, obtained through a well-defined truncation of the equations in the matter 
sector that are relevant to bound-state properties. In particular, one uses 


I{q^) = k^Q{qP, Sta") 




87r27^(l -e-9V4m?) 

ln[r + (1 + g2/A2^^)2] ’ 


(36) 


where 7m = 12/(33 — 2Nf) [typically, Nf = 4], Aqcd = 0.234 GeV; r = e^ — 1, m* = 0.5 GeV. 
A large body of observable properties of ground-state vector- and isospin-nonzero pseudoscalar 
mesons (and even numerous properties of the nucleon and A resonance) are practically insensitive 















Figure 14. The RG-invariant product d{q^) put together according Eq. (l3T]l . using the results 
of Fig. [TS) one practically obtains one single curve for all different values of the renormalization 
point ^ 



[GeV^] 


Figure 15. Comparison between the SDE-derived effective interaction obtained from the gauge- 
sector [Eq. (1351) . solid, black] with those obtained using Eq. ([36]) : pale-red band: RL truncation, 
pale-green band: DB truncation. The bands denote the domain of constant ground-state physics 
0.45 < oj < 0.6, with the curve within each band obtained with uj = 0.5 GeV. 

to variations of cj G [0.4,0.6] GeV, so long as = Doj =const. The value of ? is typically chosen 
to obtain the measured value of /j^; however its value depend on the truncation scheme used 
when deriving ([36]l . The usual rainbow-ladder (RL) truncation [83] yields ?aL = 0.87 GeV, while 
the improved truncation scheme (DB) of [76| gives ?db = 0.55 GeV. 

As can be seen in Fig. [151 the direct comparison between the T(g'^) dehned through Eq. (I3^ 
and Eq. (l36l) reveals an excellent agreement. This seems to suggest that a QCD-derived 
description of measurable hadron properties is well within our reach. 












7. Conclusions 


In this presentation we have reviewed a large body of developments related to the study of 
the off-shell Green’s functions of pure Yang-Mills theories and QCD, using SDEs and the lattice 
simulations in the Landau gauge. The general picture that surfaces may be summarized as 
follows: 

i) Gluons just as quarks acquire a dynamical mass, which is controlled by the corresponding 
“gap equation”. 

a) Ghosts remain massless, and have a finite dressing function that can be accurately obtained 
from the corresponding SDE. 

in) The quality of the ingredients entering into the quark gap equation is steadily improving, 
giving rise to phenomenologically successful quark masses. 

iv) K significant step has been taken for bridging the gap between nonperturbative continuum 
QCD predictions and hadron phenomenology. 

Even though a lot of work remains to be done in all directions, Richard Eeynman’s famous 
quote from the first of his seven “Messenger Lectures” at Cornell University in 1964 comes to 
mind: Nature uses only the longest threads to weave her patterns, so that every small piece of 
the fabric reveals the organization of the entire tapestry^'. We seem to have reached a similar 
point with QCD: even though the Green’s functions that have been studied are but a small piece 
of the QCD fabric, unequivocal signs of an underlying organization begin to emerge. 
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